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Collective behavior where a set of elements interact and generate effects that are beyond the 
reach of the individual non interacting elements, are always of great interest in physics. Quantum 
collective effects that have no classical analogue are even more intriguing. In this work we show how 
to construct collective quantum heat machines and explore their performance boosts with respect 
to regular machines. Without interactions between the machines the individual units operate in a 
stochastic, non-quantum manner. The construction of the collective machine becomes possible by 
introducing two simple quantum operations: coherence extraction and coherence injection. Together 
these operations can harvest coherence from one engine and use it to boost the performance of a 
slightly different engine. For weakly driven engines we show that the collective work output scales 
quadratically with the number of engines rather than linearly. Eventually, the boost saturates 
and the scaling becomes linear. Nevertheless, even in saturation, work is still significantly boosted 
compared to individual operation. To study the reversibility of the collective machine we introduce 
the ’entropy pollution’ measure. It is shown that there is a regime where the collective machine is 
N times more reversible while producing N times more work, compared to the individual operation 
of N units. Moreover, the collective machine can even be more reversible than the most reversible 
unit in the collective. This high level of reversibility becomes possible due to a special symbiotic 
mechanism between engine pairs. 


I. INTRODUCTION 

When particles or systems interact in a special syn¬ 
chronized way, they can collectively accomplish remark¬ 
able tasks that are far beyond the capabilities of the sum 
of individual units. An electrons Cooper pair can flow in 
a solid without resistance while unpaired electrons can¬ 
not. In biology, symbiotic relationship refers to the case 
where two different species benefit from interacting with 
each other. Similarly, in this work we study the collec¬ 
tive operation of quantum heat machines. The individual 
units are functional heat machines. The collective ma¬ 
chine that emerges by introducing a special synchronized 
interaction, is also a heat machine but with superior per¬ 
formance. Furthermore, the synchronized interaction is 
based on coherence exchange between the units and is 
therefore inherently quantum. Since this interaction af¬ 
fects only coherences in the energy basis, the machines 
do no exchange heat or work with each other. In a cer¬ 
tain case we show that the machines can be divided to 
two species that exhibit symbiotic interrelation. 

A heat engine is a device that converts heat flows into 
useful work. Work may take different forms. In the quan¬ 
tum regime it may be, for example, the amplification of 
light (lasing medium). Alternatively a quantum power 
refrigerator may use microwave radiation to take away 
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heat from a device and dispose it in some ambient envi¬ 
ronment. On top of the standard uses of heat machines, 
it has been suggested to use quantum heat machines for 
non classical tasks such as entanglement creation [T] and 
state squeezing [2]. 

With the growing technological abilities to control and 
manipulate single atoms and ions, the construction of 
microscopic quantum heat machines is within reach. A 
single atom heat engine and an NMR refrigerator have al¬ 
ready been built laia. Suggestions for realizations in sev¬ 
eral other quantum systems include quantum dots Eli, 
superconducting devices mizHi], cold bosons m, and 
optomechanical systems [HHIS]. In particular, the rapid 
advancement in fabrication technology of superconduct¬ 
ing circuits may open the road to on chip heat machines. 

The study of quantum heat machines started in m 
and has become an increasingly active area of research 
in recent years (see laiiiiiisHsi] for a very partial list). 
One of the main goals in the study of quantum heat ma¬ 
chines is to discover and understand the fundamental dif¬ 
ferences between small quantum heat machines and their 
macroscopic counterparts. 

Surprisingly, quantum heat machines show a striking 
resemblance to their classical counterparts and often ex¬ 
hibit classical features. However, recent studies have iso¬ 
lated several strictly quantum features. In [52] a proto¬ 
col was presented to extract work from coherences using 
a thermal bath without changing the energy population. 
In [53] it was shown that coherent quantum dynamics can 
lead to faster cooling rates compared to classical dynam- 
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ics. In [54] it was shown that the coherence of the work¬ 
ing fluid in the energy basis gives rise to a ’thermo¬ 
dynamic equivalence principle’ of the three main engine 
types: four-stroke, two-stroke, and continuous engine. In 
the equivalence regime these machines still have entirely 
different dynamics, yet their thermodynamics features 
are the same. This result was recently extended to the 
non-Markovian regime where the machine and the bath 
are strongly coupled [56]. For other studies of coherence 
in quantum heat engines see [1I1I57H59]. 

More relevant to the present work, is the identification 
of two distinct work extraction mechanisms [54] : a coher¬ 
ent (quantum) mechanism, and a stochastic (classical) 
mechanism. The coherent mechanism operates on coher¬ 
ences in the energy basis while the stochastic mechanism 
operates on energy populations (population inversion). 
Thus to observe quantum effects in power and heat as in 
[54] . the coherent mechanism must be more pronounced 
than the stochastic mechanism. For this to occur it is 
critical that the thermal strokes be significantly shorter 
than the thermalization time. As a result the work stroke 
(unitary evolution) starts with non zero coherences and 
the coherence work extraction becomes important. For 
a study of the role of coherence in specific machines see 
[531157]. For additional studies of coherence in quantum 
thermodynamics see [sniisiiiMiisiiinom]. 

One of the finding of the present work is that even if 
only the stochastic work extraction mechanism operates 
(i.e, the device is a stochastic machine), there is still a 
quantum resource that can be exploited for other pur¬ 
poses. This quantum resource is coherence. Although 
coherence is regenerated in each cycle in the engine by 
interacting with the classical field, it does not partici¬ 
pate in the stochastic work production process. Hence, 
the residual coherence can be used and consumed for 
other purposes without affecting the performance of the 
engine. It is shown here that coherence obtained from 
one engine can he injected to a second engine, and boost 
its power output by activating its coherent work extrac¬ 
tion mechanism. The present paper studies the possibili¬ 
ties entailed in coherence alteration processes in quantum 
heat machines. These processes not only give rise to new 
types of devices, but also have profound implications on 
reversibility and entropy balance in quantum heat ma¬ 
chines. To illustrate our findings we have chosen the 
smallest and simplest heat engine model that can exhibit 
them. Nevertheless, our findings are general and can be 
applied to other more complex devices. 

Recently, several studies considered quantum heat en¬ 
gines with a working fluid that contains multiple particles 
[65ll67] . In [68] a central spin is used to synchronize the 
other particles. These very interesting studies should be 
distinguished from the present study which focuses on 
multiple engines that operate collectively. 

Section II introduces the basic unit that will be used 
to construct a collective machine. In section III the 
over-thermalization regime in quantum heat machines is 
briefly described. Then, in section IV we introduce the 


coherence extraction (CE) and coherence injection (Cl) 
processes. In section V we quantify the amount of en¬ 
tropy reduction associated with CE. Next in section VI 
the ’entropy pollution’ irreversibility measure is intro¬ 
duced and used to show the impact of CE on non col¬ 
lective heat machines. Section VII studies collective ma¬ 
chine operation where CE and Cl are used for coherence 
sharing between the engines. After studying the perfor¬ 
mance boost gained by the collective operation we give 
our conclusion in section VIII. 


II. THE BASIC UNIT - A TWO STROKE 
THREE LEVEL QUANTUM HEAT MACHINE 
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Figure 1. (a) The three-level two-stroke engine used as a ba¬ 
sic unit for the collective machine. In the thermal stroke the 
hot bath operates on levels 1 V 3 and the cold bath on levels 
1 V 2. The unitary stroke that follows reduces the energy 
of the system and extract work, (b) The engine dynamics 
can be visualized on the Bloch sphere of level 2 V 3. The 
{ay) coordinate corresponds to coherence and the change in 
the {gz) coordinate determines the change in population in¬ 
version 5p and consequently also the change in the energy of 
the system. The top dashed arch corresponds to the pure 
states shell of the Bloch sphere. The mixed state dynamics of 
the engine takes place inside the Bloch sphere. The maximal 
population inversion is determined by the previous thermal 
stroke (red-blue arrow). An constant radius curve with an 
angle 0 (green) is formed by the unitary evolution of the work 
stroke (generated by Hw). In the over thermalization regime 
the unitary stroke generates coherence and the thermal stroke 
fully erases them. Thus, in the above setup there are no ini¬ 
tial coherences in the work stroke and the device operates as 
a stochastic machine with no quantum effects. 

The basic machine unit used in this paper is shown 
in Fig la. It is a simple two-stroke quantum engine. 
Strokes are time segments in which the engine under¬ 
goes different transformations [54] . In the thermal stroke 
levels 1 and 3 are connected to a hot bath, while lev¬ 
els 1 and 2 are connected to a cold bath. After a pe¬ 
riod that sufficiently exceeds the thermalization time, 
the populations are determined by the Gibbs factors 
P 2 /P 1 = In the work stroke 

that follows, the baths are disconnected and a unitary 
operation is applied to levels 2 and 3. Work is extracted 
during the unitary operation (e.g. the unitary can be a 
TT pulse that takes the particle to a lower energy state 
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and emits an extra photon in the process). This is a two- 
stroke quantum heat engine [39]. For more information 
about various engine types see [54]. For a more general 
introduction to quantum heat machines see reviews [69]- 
[72] and references therein. These reviews also cover other 
topics of interest in the more general and rapidly grow¬ 
ing field of quantum thermodynamics. When the bath 
are disconnected (work stroke) the Hamiltonian of the 
engine shown in Fig. 1 is 

H = Ho^ H^{t) = AEe |2) (2| + /^Eh |3) (3| (1) 

+e(i)cos[(E 3 -E 2 )t](| 2 )( 3 | + |3)(2|). (2) 

Although we will use this model for illustration, our find¬ 
ings are more general and are not limited to this specific 
three-level model Hamiltonian. The interactions with the 
baths are described by Lindblad operators that depend 
on temperature. 

The engine in Fig. 1 is an Otto engine. Unlike the 
Carnot engine, the thermal stages are isochores and not 
isotherms. This means that during the thermal stage 
the Volume’ is fixed and not the temperature. In dis¬ 
crete quantum heat machine a constant volume means 
that the energy levels are not changed in time (as they 
must change, in a Carnot machine). Since presently it 
seems that isochores are easier to implement compared 
to isotherms, Otto engines are the default choice in the 
quantum regime. The efficiency of Otto engines with 
uniform level compression in the work stroke m, and 
also for the current engine, is 77 = ^ = 1 — AEc/AE^ 
where W is the work per cycle and Qh is the heat taken 
from the hot bath. We point out that thermodynamic 
bounds on efficiency refer to steady state operation and 
not to transients or single shot non periodic evolution. 
Although the efficiency in this case does not depend on 
the temperature, it is still limited by the Carnot effi¬ 
ciency T] = 1 - AEc/AEh < 1 - Tc/Th = r]c. When 
AEc/AEh > TciTh the machine starts to work as a re¬ 
frigerator rather than as an engine. 

In [14] (see [^ for a more detailed analysis) it was 
shown that the Otto engine-refrigerator crossover point 
(sometimes called the Carnot point) takes place when the 
baths cease to generate population inversion. It is impor¬ 
tant to point out that this engine-refrigerator crossover 
is valid for “standard” quantum heat machines where 
only the baths appear in the entropy balance over one 
cycle. In this paper we introduce the coherence injec¬ 
tion process. This process changes the entropy balance 
and therefore enables the machine to perform as an en¬ 
gine even without population inversion. This is related 
to lasing without inversion m, but has a more general 
quantum thermodynamic explanation based on the no¬ 
tions of active and passive states as explained later on 
(see also Sec. HLD of [54]). 

In this study the unitary is executed by a classical field 
(e.g., a laser light, an RF pulse, or a slow magnetic field). 
This is experimentally and theoretically well motivated in 
a wide range of physical systems in atomic physics. Some 
studies replace the classical field by a quantum battery 


(work repository) [36l (56] (60] [76l [77] . This is very in¬ 
teresting, and deserves further study. Nevertheless, the 
classical field framework is highly useful in many quan¬ 
tum systems where second quantization effects are not 
important. The classical field is in some sense a source of 
coherence (see, for example, [63]). However its coherence 
is so abundant that when the classical field approxima¬ 
tion is valid there is no need to discuss the degradation 
of the classical field [63]. In short, our framework is no 
different from that used in NMR or common laser and 
atomic physics setups. 


III. BLOCH SPHERE ENGINE DYNAMICS IN 
THE OVER-THERMALIZATION REGIME 

In the over-thermalization regime the thermal baths 
are connected for periods significantly longer than the 
thermalization time, so for all practical purposes coher¬ 
ences are wiped out. Consequently, the population dy¬ 
namics in the unitary stroke can be fully described by 
stochastic means (the unitary can be replaced by a dou¬ 
bly stochastic operator [54]). Since the population dy¬ 
namics in the Markovian thermal stroke can also be fully 
described by stochastic means, the quantum machine in 
this regime operates as a fully stochastic device. Ther¬ 
modynamic observables such as work, heat and baths 
entropy generation, depend only on energy populations 
and not on energy coherences. Thus, in this regime 
there is little hope of observing quantum effects unless 
the scheme is significantly modified. The ’over’ prefix 
in the term over-thermalization [54] is due to the fact 
that beyond a certain point, large thermalization rate 
degrades the performance of the engine, as it mitigates 
the coherent work extraction mechanism. This degra¬ 
dation due to over-thermalization has been previously 
observed in specific models in [mill], but in m it was 
understood as a generic effect. In this work we show that 
even in this seemingly classical-stochastic regime of over- 
thermalization, quantum thermodynamic effects can be 
observed. 

For work extraction, the exact details of the thermal¬ 
ization mechanism are not important since throughout 
the paper we assume over-thermalization. However for 
bath entropy generation accounting we do use the weak 
coupling result to relate the heat and the bath entropy 
change via dQ = —TdSbath where T is the temperature 
of the bath. The weak coupling limit is widely applica¬ 
ble in many physical systems and constitutes the stan¬ 
dard approach in quantum open systems EH). In the over 
thermalization regime, the system is arbitrarily close to 
the thermal steady state when it ends the thermal stroke. 
The term thermal steady state refers to the state of the 
machine after it was connected to the baths (both of 
them) for a period which greatly exceeds the thermal¬ 
ization time. Stated differently the thermal steady state 
is the fixed point state of the system when it is connected 
to the baths. In general it is different from the periodic 


4 


steady state of the machine which takes into account also 
the unitary stroke. 


A. Bloch sphere representation of the engine 
dynamics 

It is highly useful to plot the dynamics of the engine 
in Fig. la in the Bloch sphere representation. In steady 
state operation, the ground state population is fixed in 
time. The unitary never changes it, and therefore to 
maintain steady state periodicity, the thermal stroke can¬ 
not change it either. Thus, the whole dynamics is en¬ 
capsulated in the two excited levels and can be most 
conveniently plotted on the Bloch sphere as shown in 
Fig. lb. The z axis corresponds to the expectation value 
of the (Jz Pauli matrix and since the interaction is 
(see 0 ) the two-level “spin” will rotate in the yz plane. 
The expectation value of 2(Tz is the population inversion 
P 3 —P 2 and the expectation values of (jy is the y coherence 
|^(P 23 — P 32 ) energy basis (there is no x coherence 

in this setup). Therefore a motion along the z axis is 
proportional to energy change, and a motion along the 
y axis is proportional to coherence change. In particular 
the 2 ; axis itself indicates zero coherence. 

In the over-thermalization regime the thermal stroke 
must end on the z axis. Unitary evolution rotates the 
Bloch vector (green curve) but cannot change its radius 
(which is equal to the half the population inversion). The 
angular spread of the unitary curve depends on the inter¬ 
action strength e, and on the stroke duration. Within the 
standard rotating wave approximation a tt pulse leads to 
a TT rotation, a 7r/2 pulse leads to a 7r/2 rotation, and 
so on. The thermal stroke generates a non unitary evo¬ 
lution (red-blue arrow). In steady state operation, this 
brings back the Bloch vector to the starting point of the 
unitary stroke and closes the engine cycle. In the over 
thermalization regime this starting point is the thermal 
steady state. 


B. Weakly driven machines 60 ^ tt 

In this paper we shall assume that the units are weakly 
driven so that 0 ^ SO ir. Nonetheless, SO is finite and 
not infinitesimal. Weakly driven machine don’t exhaust 
their population inversion in each cycle and therefore, in 
some sense, don’t reach their full potential. Nevertheless, 
this does not necessarily means low power since smaller 
changes in the engine can be completed in a shorter time. 

Clearly by increasing the external driving field, more 
work can be extracted. However, in this work we shall 
assume that small SO engines are given and that SO ean- 
not he inereased beyond the given value. The question is 
what can be done to increase the work output without 
adjusting the drive. For a single unit this is not possible, 
but we show that for multiple units that operate collec¬ 
tively it is possible. In practice, ^ 7r/3 is small enough 


to start observing the boost of the collective machine. 

Weak drive may result from practical reasons. One 
reason may be the lack of strong enough coherent laser 
source in the needed frequency. Alternatively the source 
can be strong but off the atomic resonance. The large 
detuning reduced the amplitude of the Rabi oscillation 
(directly related to SO). Another motivation to use weak 
drive is to avoid excitation of other undesired mechanism 
and effects such as nonlinear multiphoton processes. 


IV. COHERENCE EXTRACTION AND 
COHERENCE INJECTION 

With the exception of exact permutation, any unitary 
that operates on a diagonal state generates coherences in 
the same basis. This is especially true for unitaries that 
are close to the identity where this is the leading effect 
generated by the unitary. For diagonal states and a small 
time step dt, the change in population is proportional to 
dt^ while the change in coherence is proportional to dt. 
However, as explained earlier, in the over thermalization 
regime this coherence is erased by the baths in the next 
stroke. 

The idea that facilitates our main results is simple. In¬ 
stead of letting the bath erase the coherence we wish to 
extract and store it elsewhere for further use. We impose 
the condition that the extraction does not change the en¬ 
ergy population as this may lead to heat or work energy 
exchanges that significantly modify the thermodynamic 
scheme. Therefore, the coherence extraction has no en¬ 
ergetic role only an entropic one. A simple way to imple¬ 
ment this extraction is to swap the engine particle with 
a particle that has the same energy levels and the same 
populations but without any coherences. We call this 
particle coherence acceptor. This coherence extraction 
scheme is illustrated in Fig 2a. Since the engine operates 
in the stoehastie regime (over thermalization), only pop¬ 
ulations matter, and this swap operation will have zero 
impact on the operation of the engine. All energy cur¬ 
rents (heat and work) will remain as they were before the 
CE. In the weak coupling limit to the baths (that leads 
to the Markovian dynamics), the bath entropy genera¬ 
tion is determined by the heat via ASbath = ~Qt^ {T is 
the temperature of the bath). Therefore ASbath is also 
indifferent to the presences or lack of coherences at the 
beginning of the thermal stroke. 

To graphically describe CE two Bloch spheres are 
needed: one for the engine and one for the acceptor par¬ 
ticle. As shown in Eig. 2b, at U, just before the thermal 
stroke, a swap operation starts between the engine and 
the acceptor (golden arrows). When it ends at ^ 2 , the 
acceptor has coherence while the engine reaches the ^ 
axis and looses all coherence. Since both particles have 
the same populations, the population inversion {z coor¬ 
dinate) is fixed as well during the swap. As a result, the 
system energy does not change and the CE operation 
does not involve any heat or work exchange between the 
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Engine Coherence acceptor 


Figure 2. Introducing the coherence extraction (CE) process, 
(a) A density matrix illustration. The thermal stroke (bot¬ 
tom left of (a)) sets the engine in a diagonal state. Then, 
the unitary stroke generates coherences and some population 
changes in level 2 & 3. Next, the engine particle is swapped 
with an acceptor particle that has the same energy levels and 
the same population. This swap does not involve heat or work 
exchange. The engine state becomes diagonal but the engine 
performance are not affected by this. On the other hand, the 
gained coherence of the acceptor can be used to boost the 
performance of a different engine that operates with slightly 
different temperatures. This shows that even in the stochastic 
regime quantum heat engine have quantum byproducts that 
can be used, b) Illustration of CE in Bloch sphere represen¬ 
tation (see text). 

system and the auxiliary particle. 

Later on coherence extraction will be used together 
with its complimentary process, coherence injection, to 
generate interaction between machines that lead to col¬ 
lective operation. However, before using CE for con¬ 
structing a new type of quantum heat machines, it is 
worthwhile to discuss some immediate implications of CE 


Benefits of coherence extraction 
A. Cooling via CE 

Coherence extraction adds coherence to some external 
particle without affecting the operation of the engine. 
What is this coherence good for? First, in terms of en¬ 
tropy, adding coherences while keeping the diagonals the 
same, means that the von Neumann entropy of the parti¬ 
cle has been reduced (the state becomes purer). This is a 
form of cooling but without changing the average energy 
of the particle. Since the entropy of the auxiliary parti¬ 
cle has been reduced and the entropy production in the 
baths remained the same, the total entropy production 
in the world (see entropy pollution later on) has been 
reduced. After we quantify it in Sec. VI it will become 


clear that this can be a dramatic effect with important 
consequences. 


B. State activation via CE 

In terms of energy, the added coherence is also impor¬ 
tant since the state becomes more active (less passive) 
than it was before. Passive states are diagonal states in 
the energy basis whose populations are decreasing with 
energy [79H8T] . Stated differently, passive states have 
no coherences and no population inversion in the energy 
basis. All other states are termed ’active states’. The 
reason for the names ’active’ and ’passive’ is the follow¬ 
ing. It is impossible to extract energy from a passive 
state by a unitary transformation if the Hamiltonian is 
returned to its initial value after the transformation is 
completed (e.g. a tt pulse). On the other hand, some 
energy can always be extracted from active states using 
unitary operation until they become passive 179]. 

Even if the state had some population inversion so that 
it was not passive to begin with, adding coherences will 
make it more active. To see this, consider a two-step 
process where in the first step permutations are used to 
extract the energy from the population inversion. This 
stage extracts work which is equal to the work extracted 
from the same original state but with no coherences. 
However, since permutations do not alter the magnitude 
of coherences, the state remains active. In the second 
step, a rotation is performed to bring the state to a di¬ 
agonal form (whose energy is guaranteed to be lower). 
Thus the coherence extraction make the auxiliary state 
more active and therefore increases the energy that can 
be extracted from it. In some cases the acceptors need to 
have some population inversion, that is, they need to be 
active to begin with. However, one can show that with 
help of two TT pulses it is possible to use a passive ac¬ 
ceptor in this case as well (without an additional energy 
cost). 

At this point it is clear, at least qualitatively, that CE 
can cool and activate a state simultaneously. We point 
out that in principle, in other processes, it is possible to 
cool without activation and to activate without cooling. 
Here, however, these two good things happen simulta¬ 
neously. Naturally the question that comes to mind is 
whether these gifts come for free, or are there some ad¬ 
ditional consumed resources? As explained earlier CE in 
over-thermalization does not change the heat and work 
flows of the engine. Nevertheless, there could be other 
non energetic resources. The immediate suspects for ex¬ 
tra resources are the acceptor particles, since they have 
to be prepared in a specific state. 

In [82] we have presented a new way of keeping track 
of resources in a very general setup, and showed that the 
acceptor particles in CE are not resources but separate 
objects that perform a different task from that of the 
engine. This extra task uses the coherence byproduct 
of the engine as fuel. Hence the cooling and activation 
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are not exactly free, but are carried out without extra 
additional resources on top of what the engine is already 
consuming (heat). The curios reader is encouraged to 
follow the logic of resource accounting in Sec. VIII of 
[82] (vl). However this accounting is not needed for the 
collective machine analysis studied here. 

Before doing the entropy balance for the CE process 
let us introduce another useful process and then do the 
entropy balance for both processes together. 


COHERENCE INJECTION 

The coherence injection (Cl) process shown in Fig. 3a 
is the time-reversal, complementary process of coherence 
extraction. A ’coherence donor’ particle with coherence 
is swapped with the engine particle that has no coher¬ 
ences. The donor particle and the engine particle have 
the same energy levels and the same energy populations. 
Thus, like CE, Cl does not involve any energy exchange. 
However, in contrast to coherence extraction this process 
is carried out before the unitary stroke. The injected co¬ 
herence allows to extract more work for the same uni¬ 
tary operation. This is easily understood in the Bloch 
picture. Consider a unitary rotation with a small angle 
60 TT. If there is population inversion but no coher¬ 
ence, the change in the projection on the 2 ; axis (work) 
is 1 — cos 60 = ^{60)‘^ -1- O[{60)^]. On the other hand, 
with coherence there will be a linear term as well. If 
the maximal amount of coherence is injected [83] before 
the rotation, then the Bloch vector starts on the equator. 
This time the small rotation will lead to an energy change 
proportional to 60. Alternatively stated, the injected co¬ 
herence enables to get larger change in the populations 
(z coordinate) for the same rotation 60. 

Another use of coherence injection is to change the 
character of the machine without changing the tempera¬ 
tures. As mentioned earlier without Cl the engine in Fig. 
1 must satisfy AEc/AEh > TcjTh in order to have pop¬ 
ulation inversion and extract work. In Fig. 3b we show 
that if coherence is injected, it is possible to extract work 
(have an engine) even when there is no initial population 
inversion. Once the Bloch vector is not pointing to the 
south pole, i.e. it is an active state, it is possible to ap¬ 
ply a rotation that will decrease the energy of the system 
{z coordinate) and extract work by rotating the vector 
towards the south pole. Note that the efficiency is still 
given by 1 — AEdAE^. In this particular engine Cl al¬ 
lows do more work but this is compensated by taking 
more heat from the hot bath so the efficiency stays the 
same. As pointed out in [ZSIIHI], coherence can be used 
to extract energy from a single bath. This is fully con¬ 
sistent with the second law when taking into account the 
coherence that has to be supplied to the machine. 

On its own coherence injection is interesting but not 
very practical. In Cl it is clear that coherence donors 
are a consumed resource that has to be accounted for. 
Thus, in the present paper Cl is used only as part of a 



Figure 3. (a) Coherence injection (Cl). Cl is the time inver¬ 
sion of CE. Potentially it can dramatically boost performance 
as the unitary stroke can produce more population change if 
it starts with some initial coherence. Due to the complexity 
of preparing the coherence donor particles, this process will 
only be used as a part in a larger device where the donors are 
generated automatically by other parts in the device. In fact, 
one engine will be used as a donor for a different engine that 
will be used as an acceptor. Cl can significantly extend the 
working regime of engines, (b) Shows an engine operating 
without population inversion. This engine can operate also 
when Tc = Th (a single bath) since the influx of coherence 
prevents a conflict with the second law. 

larger deviee. In this larger deviee the donors are ereated 
as a byproduet of a different element in the deviee. This 
scheme is very different from other proposals that use 
“coherent baths” to supply the machine with coherent 
(non-thermal) fuel (“phasonium” |84l|85| and “squeezed 
baths” HilHZ]). In [53| it is assumed that coherence is 
somehow added to the setup. Coherence injection is a 
way to add this coherence without energy costs. 

One may suspect that in such a closed device that has 
only work and thermal baths terminals, all benefits of 
Cl and CE will vanish and the device will operate as a 
standard thermal machine. However we show that this 
is not the case and that the internal Cl and CE lead to 
observable quantum thermodynamic signatures, and to 
significant performance boost. 

V. THE SECOND LAW WITH ACCEPTORS 
AND DONORS 

Earlier it was qualitatively shown that the engine re¬ 
duces the entropy of the acceptor particles (one particle 
per cycle). In this section we want to quantify this effect, 
and show that the scaling properties of this cooling effect 
has a very significant consequence to the theory of heat 
engines and their reversibility. For this purpose we find 
the following two identities extremely useful: 

S{p 2 ) - S{pi) = -tr[{p 2 - pi) \npi] - D{p 2 \\pi), (3) 
S{p2) - S{pi) = -tr[{p2 - pi)\np2]+ D{pi\\p2), (4) 

Where S{p) is the von Neumann entropy and D{p 2 1 |pi) = 
tr[p 2 (lnp 2 — Inpi)] >0 is the quantum relative entropy 
[55] . Using @ for the acceptors, and @ for the donors 
we readily obtain 

^Sacpt = -D{plilT\\Pacpt) = -C{p^), (5) 

^Sdnr = D{pdnr\\p"dnrl = +C{pdnr), (6) 
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where is the density matrix after the unitary stage and 
Pacpt 5 pdnr are the state of the acceptor and the donor be¬ 
fore the CE and Cl take place. The coherence measure C 
is given by C{p) = D{p\\diag{p)) > 0. C was identified in 
[89] as a proper coherence monotone. Furthermore, this 
form appeared in m, in the context of thermodynamic 
transformations. 

Our first main finding is that the total entropy AStot = 
ASbaths + ASacpt + ASdnr after one cycle of the machine 
is 

ASfot ~ ^{Pw\\Peq) t ( 7 ) 

where are the populations after unitary stage and Peq 
stands for the thermal steady state population. The en¬ 
gine entropy does not appear in the entropy balance be¬ 
cause in steady state the engine state is periodic so after 
one cycle AS engine = 0 - Since the relative entropy is al¬ 
ways positive (or zero), equation 0 contains within the 
2nd law. However, Q is stronger than the second law 
AStot ^ 0; the exact equality can be used to study quan¬ 
titatively the entropy generation and its scaling proper¬ 
ties. It is interesting that after including the quantum 
effects of CE and Cl, the second law can be accurately 
described by populations only. To prove ^ we write: 

ASacpt _|_ /^^^gdnr ^ 

[S(^Pw) S{Pacpt)] T [^{Peq) ^{Pdnr)] ~ 

S{peq) - S{pacpt) = S{peq) “ S{p^), ( 8 ) 

where we have used the fact that pdnr is related to p^ via 
unitary transformation and therefore S{pdnr) = S{pw)- 
Using (Isl) and setting pi = diag{pw)^ p 2 = diag{peq) in 

0: 

qacpt I ^i^dnr _ 

^ ^X P&qd Pw^j^^^PeqP 'S D{pw\\Peq) 

3 

~ ^ ^X Peqp ~ Pw,j){~AEj/Tj + Z) + D[pyj\\peq) = 

3 

Qh/Th E Qc/Tc E D{pw\\Peq)^ ( 9 ) 

where Z is a normalization factor that comes from peq. 
Since the first two terms on the right hand side of ^ are 
equal to -ASbaths, 0 immediately leads to Q . 

To appreciate the dramatic implication of ([7), we study 
entropy generation in engines where the unitary opera¬ 
tion only slightly changes the populations (very small 
rotation) with respect to the thermal steady state. Let 
us denote the population change as Sp = p^j — Peq • the 
expression for the baths entropy is: 

^Sbaths = ^ lnpe9,i- (10) 

For small Sp the expansion of the relative entropy is: 

AStot = D{Pw\\Peq) = 'Zs—~ ~ 

^^2peq,j Qpigj 

( 11 ) 


ASbaths contains a linear order in 6 p but as we see 
from (11), AStot is only quadratic in 6 p. This means 
that in coherence extraction, a remarkable cancellation 
of the entropy generation leading order 0{5p) ( [1Q| ) takes 
place. Thus, the cooling effect of the acceptors is not a 
minor effect even with respect other entropy modifying 
elements (baths). It is a leading order effect and as we 
shall see in the next sections, it can lead to new types 
of machines and to new reversibility limits. Another im¬ 
portant property of 0 is that it holds regardless if there 
is a donor or not (acceptor is necessary). However, the 
magnitude of Sp is affected by the donors coherence as 
explained earlier. 


VI. ENTROPY POLLUTION AND A NEW 
REVERSIBLE LIMIT FOR OTTO ENGINES 


Entropy generation gives some information about re¬ 
versibility. However, to a certain extent, only the limit of 
zero entropy generation is perfectly clear. In particular, 
entropy generation alone is not sufficient for comparing 
the reversibility of two engines. For example, one en¬ 
gine has AStot = 1 and work W = 1 and the other has 
AStot = 1 and W = 1000. Clearly, the second engine is 
more reversible. To produce the W = 1000 as the second 
engine a thousand engines (or thousand more cycles) of 
the first kind are needed. The total entropy production 
of a thousand engines of the first kind is AStot = 1000, 
which is clearly less reversible compared to the second 
engine with AStot = 1 for the same amount of work. 
This motivates us to define the entropy pollution (EP) 
irreversibility measure as: 


EP = 


AStot 

Wtot ’ 


( 12 ) 


where the quantities are calculated for a whole cycle in 
steady state. Note that in m the same measure was 
used to quantify the charging quality of a battery. Here, 
however, we use (12) to quantify the reversibility of the 
whole setup including all the elements. One might be 
tempted to suggest r^c — ^ as a measure of irreversibility, 
but the problem is that when there are multiple baths 
as in the collective machine the Carnot efficiency is not 
uniquely defined. Depending on the bath connectivity 
there could be many different Carnot machines with dif¬ 
ferent efficiencies. In contrast, the definition (12) is free 
from this ambiguity. 

For classical or quantum engines without coherence in¬ 
jection or extraction: 


EP, 




-iW-Qh)lTo + Qhm 

W 

(1 - 7i)/To - 1/Th 


v 


1 Pc-P 
Tc p 


(13) 


For the Carnot efficiency the entropy pollution ( |13| ) yields 
zero as expected. This general expression for 












shows that as long as the efficiency and temperatures 
are kept the same, the entropy pollution is the same. 
It does not matter if there are multiple units operating 
independently in parallel or if the value of Sp is changed. 
Just as an example, for the Otto engine in Fig. 1 77 = 1 — 
so we verify that without CE the entropy pollution 
does not depend on Sp (i.e. 0(Sp^)). Next we study the 
case of an engine with CE. Using we find 


EPc/C£; = 


^{Pw I \Peq) 

EUejSpj 


= 0{6p). 


(14) 


where we used W = ^j=i EjSpj = 0{Sp) and (11). This 
result has a profound implication as we illustime next. 
Let our goal be the production of Wq work using one en¬ 
gine. This Wq is associate with some population change 
per cycle Spo. Now instead of doing Spo per cycle we 
only do Spo/N per cycle and perform N time more cy¬ 
cles compared to the previous scenario. Clearly the same 
amount of work Wq is obtained in both scenarios. Yet, 
in the smaller steps case, the engine operates closer to 
the thermal steady state Peq. Although this is similar in 
spirit to a quasi-static process, it is easy to verify from 
(13) that without CE, the entropy pollution is not affeeted 


by this splitting to small steps. Even though the system 
is now always closer to the thermal steady state, the ir¬ 
reversibility (entropy pollution) has not decreased. In 
stark contrast, according to ( p^ with CE we get 


ITo 


(15) 


By taking the quasi-static limit N ^ oc we now obtain 
that EFce 0- This is a reversible limit as no entropy 
generated in the total setup (acceptors-Fbaths+engine). 
At first, this result seems strange. This Otto engine pro¬ 
duces finite work and its efficiency is less than Carnot 
(77 = l—EcjEh)^ and yet the entropy pollution is zero, i.e. 
it is reversible. This seems odd as the Carnot’s theorem 
states that two-bath reversible machines must operate at 
Carnot efficiency. Indeed, without CE the second law 
permits only one type of reversible machine - the Carnot 
machine. However, the inelusion of CE ehanges the en¬ 
tropy halanee and enables the existenee of other types of 
reversible maehines whose effieieneies are lower than that 
of Carnot. In fact, when expanding the acceptors entropy 
for small Sp one finds that the bound predicted by the 
second law AStot = ASbath + ASacpt = 0 leads in the 
limit J 79 ^ 0 to the observed efficiency 1 — AEdAE^ 
and not to the Carnot efficiency. Once again we have full 
consistency with the second law. 

Note that this reversible limit does not have to be slow. 
Instead of doing N times more cycles, it is possible to 
have N engines working in parallel for the original num¬ 
ber of cycles. Other aspects of reversibility will be dis¬ 
cussed in a different publication. 

In the Sp ^ 0 limit entropy is not generated. It is sim¬ 
ply redistributed between the baths and the acceptors. In 



a) 


Sp„ / N. 

/ 
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I 


N 
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rr T 


TOT 


w 

rr ^ 


TOT 


> No improvement in 
entropy pollution 



With CE; 


AS 


TOT 


w 
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N 


■^0 


^ CE induced 
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Figure 4. (Top) Splitting the work done in one cycle Wtot 
into N many cycles, makes the system stay closer to the ther¬ 
mal steady state. However, as shown in (13) the entropy 
pollution irreversibility measure is not reduced or affected by 
this splitting (see (a)). In contrast, when CE is added to the 
scheme (see (b)), the entropy pollution goes to zero in the 
limit of large N. This enables to achieve reversibility even 
though the efficiency is lower than the Carnot efficiency. As 
explained in the text this is fully consistent with the second 
law. 


the scenario studied above the heat and work flows were 
not affected by the coherence extraction. To avoid un¬ 
necessary entropy generation coherence was stored in the 
acceptors that act as coherence repositories. Although 
this is interesting from an informational point of view, 
our main goal in this paper is to show new ways of har¬ 
nessing coherence to boost performance. Our approach 
to achieve this performance boost is to construct a col¬ 
lective operation scheme as discussed in the next section. 


VII. COLLECTIVE COHERENT ENGINES 

As we have shown, CE in quantum heat engines re¬ 
duces the entropy production but it does not affect the 
power performance of the engine. It enables other tasks 
like cooling and activation for free, but the main function 
of the engine, power production, does not benefit from 
the free coherence production. Now we present a scheme 
where CE and Cl are used to simultaneously boost power 
(not efficiency) and reversibility. This scheme has an¬ 
other very appealing property; it does not use any exter¬ 
nal aeeeptor and/or donor partieles as before. It is based 
solely on coherence sharing between engines. When two 
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engines interact one engine perform as donor and the 
other as an acceptor. 

Earlier we introduced the coherence injection process. 
Despite the potential power boost it offers, the prepara¬ 
tion of a stream of coherent donor particles makes the 
scheme less appealing. Fortunately, as shown earlier, en¬ 
gines produce coherence that can be harvested and used 
without any energetic cost. It is natural, then, to use 
one engine as a donor and pass its coherence to a second 
engine that will serve as an acceptor. The first engine 
experiences coherence extraction and the second engine 
experiences coherence injection. The second engine gets 
a power boost by starting the unitary with some coher¬ 
ences as explained earlier and in [54]. The total power 
will be the power obtained from the two engines (the 
first stage also produce power). This process can be re¬ 
peated for multiple stages until the is no more coherence 
to harvest from the last stage. In what follows we study 
two different regimes and compare to individual opera¬ 
tion without coherence sharing. 


A. Quadratic power boost for collective machines 

Our collective coherent engine scheme is shown in Fig. 
4a. The first engine donates its coherence to the second 
stage engine. The temperatures of the second engine 
baths are chosen so that resulting thermal steady con¬ 
stitutes an acceptor for the first engine. It is a perfect 
synergy where one stage prepares what the other stage 
needs while all units produce work at the same time. 
In addition, later on in the cycle the second stage en¬ 
gine donates its coherence to a third-stage engine, and so 
on. Each donor-acceptor interaction between the engines 
is carried out by a full swap of the two engine parti¬ 
cles. Thus only engine particles are participating in this 
scheme and there are no external donor-acceptor parti¬ 
cles. The question of Cl extra resources is irrelevant in 
this device. 

Due to this process of coherence transfer, each engine 
continues the unitary evolution where the previous stage 
has stopped (green line in Fig 4a and 4b). Thus, instead 
of many copies of 50 rotations as happens without CE- 
CI, the whole system performs 50 unitary as shown in Fig 
4b. It is like an assembly line where eaeh engine does its 
part. We denote by O = NdO the collective angle. For 
initial population inversion dpo the total change in the 
population inversion is 5p = dpo{l — cos Q) = 2dpo sin^ ^ 
and therefore the work of the collective device is 

Wcoii. = (Eh - Ec)dpo sin^ (16) 

Next, we wish to compare the collective operation to the 
operation of N independent units with the same 50. Now 
that there is no CE and Cl there is no need for using dif¬ 
ferent temperatures for different engines. Nevertheless, 
for the fairness of comparison, we do not allow the stan¬ 
dalone units to be connected to any cold and hot tem¬ 
peratures. The temperatures must be chosen from those 


available to the collective machine. This principle will be 
used in other comparisons later on. 

We choose the unit that produces the maximal amount 
of work when working in a standalone mode, and set the 
other N — 1 units to work with the same temperatures (N 
independent copies of the same engine-bath setup). The 
engine that produces the maximal amount of work as a 
standalone unit is the first one, because it has the highest 
initial population inversion. Thus the total power in the 
standalone work optimized (SWO) case is 


Wn,swo = N{Eh - Ec)dpo sin ^ 


n/N 

2 


The work boost for Q < tt: 

WcoiL ^ ^ ^ ^ 

Wn,swo N SW? 



(17) 


(18) 


When Q ^ TT the collective work boost is simply N (the 
coefficient is unity and not 4/7r^). Alternatively stated, 
while in the standalone ease the work for a fixed dO seales 
like the number of units N, in the eolleetive ease the work 
seales like N‘^. 

We point out that in the performance comparison stud¬ 
ied in this paper the ’unitary strength’, 50^ of each units 
is fixed. We do not compare machines with different val¬ 
ues of 50. Instead, we ask the following question: given 
that N machine with fixed 50 are available, is it pos¬ 
sible to boost their performance by making them work 
collectively? 


B. Saturation and boosted linear scaling of 
collective machines 


The expression (16) for the collective work is bounded 
and reach its peak value at Q = tt. On the other hand, 
the expression for the non interaction SWO units always 
scales like N. Does this means that for some large enough 
value of N the collective machine will produce less power 
than N SWO machine? not quite. Despite the satu¬ 
ration, collective operation still give rise to significant 
power boost. 

Let M = it/ 50 be the number of units needed for 
a Q = TT collective rotation. We have at our disposal 
N ^ M units, each of strength 50 (which is fixed and N 
independent). There is no point in letting more than M 
units interact together. For example 2M units operating 
collectively will generate a 27r rotation that amounts to 
zero work. Hence, we divide the units into sets of M 
units (one set may contain less than M units). Each set 
is operated as a collective machine. For the complete 
sets it holds that Q = tt, so according to ( [T^ their power 
amounts to {Eh — Ec)dpo [A/MJ where := Eloor{-). 
Adding the power from the single incomplete set we get 


Wcoll,N>M = 

{Eh - Ec)dpo{[N/M\ + sin^(i?em(A, M)50/2)} 

( 19 ) 
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Figure 5. (left) In our collective heat machine scheme coher¬ 
ence is extracted from one machine and injected to the other. 
The resulting “collective unitary” of angle Q produces more 
than N times more work compared to the standalone opera¬ 
tion (right). The work boost is N even when comparing to N 
standalone units that produce the maximal amount of work 
(SWO units). This means that for a fixed 60 the collective 
work scales quadratically with the number of units rather than 
linearly as in the standalone operation. 


where Rem is the remainder. For large N and small 60 
{60 7r/5 or M 5): 


Wcoii^N^M ^ N/M _ 4M 
Wi^SWO 60^14: TT^ 


( 20 ) 


Comparing this result to the SWO scaling = N 

we observe that there is a ^ power boost even though 
the behavior is no longer quadratic for N > M. This is 
illustrated in Fig. 6 for M = 20. The graph aim to show 
the N > M behavior. The red dot show the actual power 
of the collective machine as give by expression (19). Ini¬ 
tially the power grows quadratically (blue curveJT^ After 
that the power follows the boosted linear scaling (green 
curve) given by (20). For comparison, the orange curve 
shows the linear power growth of non interacting units. 
For smaller 60 the quadratic scaling regime is larger. 



Figure 6. When the number of units N exceeds the number 
of units needed to complete a tt rotation in the Bloch sphere 
M, the boost changes from quadratic (blue) to linear (green). 
Nevertheless the linear scaling of the collective machine is 
still boosted by a factor 4M/7r^ (for M >~ 5) compared to 
the standard linear scaling N (orange) of the non interacting 
SWO units. Depending on 60, the quadratic regime can be 
larger. 


are not in conflict. Both of them increase with the num¬ 
ber of units. This is due to a unique symbiotic mechanism 
that takes place in the Q = tt collective machine. 

In the O ^ TT case all units are close to the north pole, 
so the entropy pollution of the collective machine is very 
similar to that of a standalone machine with temperature 
that corresponds to maximal population inversion (north 
pole of the Bloch vector). 

In addition, in the Q < it collective machines there is 
some residual unused coherence since the Bloch vector 
does not return to the 2 : axis at the end of the collective 
unitary {ft rotation). Since there is no external acceptor 
to store this coherence it is being erased by the baths. 
Consequently, the bath entropy generation is larger com¬ 
pared to the case where coherence is stored. 

For a fair performance evaluation of the Q = ir engine, 
it is important to take the right comparison reference. We 
shall make two different comparisons. First we compare 
the O = TT case to the standalone work optimized case 
studied in the previous section. Next we shall remove the 
work optimization constraint and compare to the most 
reversible unit in the collective. First we start with the 
scaling of the entropy pollution in the Q = ir case. 


C. Extraordinary reversibility via symbiotic 
relations at Q = tt collective angle 


D. Entropy pollution of the Q = tt machine 


In this section we study the intriguing case of N60 = tt. 
In terms of work boost we simply need to set = tt in 
(18). However, in this section we are interested in the en¬ 


tropy pollution of the collective device. Surprisingly, we 
find that the collective machine can be more reversible 
than the most reversible unit in the collective. Moreover, 
in this machine work yield boost and reversibility boost 


In our collective scheme all units but the last one, have 
acceptors (the next stage engine) that remove the re¬ 
maining coherence after each engine complete its unitary 
stroke. Hence, in the general case, the entropy genera¬ 
tion in all units but the last can be evaluated using 0 . 
However, when Q = ir the last machine ends at the south 
pole with no coherences so there is nothing to extract. 
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Figure 7. The Q = tt collective machine shows a remarkable 
symbiotic interaction between units in the upper and lower 
hemisphere. The coherence the lower units require in order 
to operate as engines (they have no population inversion) is 
supplied by the upper units. In return the lower units produce 
negative bath entropy that in the limit of large N compensates 
for the positive entropy generation of the upper units. The 
entropy generation in the baths of each units is shown in the 
rightmost graph. In the limit of large N the graph becomes 
antisymmetric and the total entropy production goes to zero. 
In this limit the Q = tt is reversible. 


Thus, in the Q = tt case we can safely use 0 for all 
engines including the last one, and get that the total en¬ 
tropy production is 


N 


AStot = ASZlL = 


( 21 ) 


i=l 


Peq = e , e }/(l + e 


+ e 


AEh 


( 22 ) 


where we used = P^w • This p^eq^^ = P^w con¬ 

dition determines the temperatures of the i + 1 baths 
- it is the basic population matching in CE and Cl 
that assures that there are no energy exchanges be¬ 
tween the units. In addition, AS'tot = ASbaths fol¬ 
lows from the fact that in this system there are no ex¬ 
ternal particles so in a periodic process the total en¬ 
tropy of the elements can only be modified in the baths. 
For, small population change 6 p associated with small 


60, we can write: D{p, 


(*+i) 

'eq 


I 


_ 1 

) 2 ^7 = 




PeqU) 


and 


Spi = dpo sinOiSO, where j is a level index and i is an en¬ 
gine index. Using dO = ir/N we have ASbaths = AS tot = 

1 '^3 

2-^i=l 2 A/j = 


piqij) is a function of di we get ASbaths = AStot = 
F{0i)60. The sum can be approximated by an in¬ 
tegral for large N. Since the integral does not depend on 
N we get that AS^F^ cc 1/N. The total work in this case 
is dpo{Eh — Ec) regardless of N, so finally, the entropy 
pollution scaling is: 


Spj 

PeHj) 


_ 1 'spS dpQsin OiSO^ qfnpp 

- 2 2.,=2 • ^ince 


(23) 


Comparison to the standalone work optimized ease 


We saw that if the N units that constitute the collec¬ 
tive engine, work independently and the temperatures are 
chosen to maximize the work, their total work is N times 
weaker compared to the collective machine (see @). 
One might suspect that this enhancement is accompa¬ 
nied by degradation of reversibility, i.e. in an increased 
entropy pollution. In the standalone work optimized case 
the entropy pollution is the same for all units. It is given 
by (13) with the temperature and Carnot efficiency that 


fits the north pole of the Bloch sphere. The exact ex¬ 
pression is not important. What matters is that EPgwo 


does not depend on the number of units N. Using (23) 
we immediately get 


EPn= 


0=7r 




0{^/N) 

0{NO) 


= 0{1/N) 


The conclusion is that the Cl = tt collective machine pro¬ 
duce N times more work, and at the same time it is N 
times more reversible compared to the standalone work 
optimized case. In addition this shows that even though 
there are no external auxiliary particles, CE and Cl are 
still highly useful as they can increase work yield and re¬ 
versibility. 

Note that the comparison above is very strict. If the 
coherence exchange is not performed but the tempera¬ 
ture of each unit is the same as in the collective machine, 
then all the lower hemisphere units will turn into refriger¬ 
ators and consume work instead of producing work. As a 
results, the total work output in this ’no Cl CE’ case be¬ 
comes zero for large N (the cancellation is not perfect for 
small values of N). Of course in such a case there is little 
motivation to use units that do not produce work. This 
is why we chose to compare it to the power-maximized 
standalone case. Next, we compare it to the standalone 
EP optimized case. 


Comparison to the most reversible unit 

If we look on the individual units operating at differ¬ 
ent temperatures needed for the collective machine, the 
unit with the lowest entropy pollution is unit number 
\N/2] := Ceiling(A/'/ 2 ), the one that starts just above 
the equator m- This unit has the lowest entropy pol¬ 
lution since for large N its Carnot efficiency approaches 
the actual efficiency 77 c ([A/2]) rj = 1 — AEdAE^ (see 
Now we wish to compare the collective machine 
to the most reversible unit (in standalone operation) in 
the collective. Following the same logic as in the SWO 
case, N copies of this most reversible unit are chosen as a 
comparison reference. We denote this reference by SEPO 
for ’standalone EP optimized’. We once again point out 
that for the fairness of comparison, we are only allowed 
to use bath temperatures that are available to the col¬ 
lective machine. It is simple to show that this time the 


EPn=n oc 1/N. 
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ratio between the work of the collective Q = it machine 
and the SEPO scales like N‘^. However, when compar¬ 
ing the entropy pollution of the two scenarios one finds 
something surprising. 

After some algebra one can show that the entropy pol¬ 
lution of the SEPO scales like 1/A^, just like the ft = tt 
collective entropy pollution (23). In Eig. 8 the ratio 
EPcoZ/./EPsepo is plotted for different minimal temper¬ 
atures. As expected from the scaling arguments above, 
the ratio saturates for large N. The most striking fea¬ 
ture in this figure is that the collective machine can be 
more reversible (less polluting) than the most reversible 
unit (the plotted ratio is smaller than one). The second 
feature that stands out in Eig. 8 is the dependence on 
the parity of units number. This parity dependence is 
not coming from EFcoii. but from the SEPO reference. 
Eor even N the most reversible machine is at an angle ^ 
from the equator while for odd N the angle is Eor 
large N this means a factor of two in the rjc — V factor 
that appears in (13). The figure shows that when the 


cold temperature is comparable (or higher as we see nu¬ 
merically) to the cold gap, the collective machine is more 
reversible than the SEPO regardless of parity. When the 
low temperature is below the cold gap we observe that 
only the even N values lead to superior reversibility of 
the collective machine. 

Of course this small improvement shown in Eig. 8 is 
not a performance boost. The performance boost is in 
the power as explained in the beginning of this section. 
Nevertheless, it is remarkable that the collective machine 
can be more reversible than the most reversible unit - 
even by a little. 

The two comparisons shown in the upper table in Eig. 
7, can be combined by considering the figure of merit 
W/EF that can be understood as work multiplied by re¬ 
versibility measures 1/EP. Repeating the same argument 
as before it is easy to verify that the W/EF measure of 
the ft = TT machine is N‘^ better compared to both SWO 
and SEPO standalone operation. 

In what follows we present a different point of view 
that unravels the interesting dynamics that leads to high 
reversibility and high work yield in the ft = ir collective 
machine. The rightmost graph in Eig. 7 shows the bath 

entropy generation ASbaths,i = -Qh^- Qc'’ 
of each engine (red circles N = 6, blue squares N = 
20). In the lower hemisphere there is a negative entropy 
production. Clearly, from this becomes possible only 
if there is CL The source of the injected coherence are 
the engines in the upper hemisphere. As the number of 
engines increases the graph becomes more antisymmetric 
and the total sum of entropy generation in the various 
baths becomes closer to zero. The reason for the almost 
exact, leading order, cancellation of the i baths entropy 
generation with that of N -\-l — i baths (antisymmetr;^, 
can be understood in the following way. Equation (§-0 
lead to 


^^baths,i — Ci^Pdnr^i) C{,Pw,i) f^ijPw ,i\\Peq,i) • 
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Figure 8. Comparing the entropy pollution of the collective 
machine to its most reversible element in standalone units 
(SEPO). The graph shows that the reversible eolleetive ma- 
ehine ean he more reversible than the most reversible unit that 
eomposes it (pollution ratio smaller than one). The parity de¬ 
pendence is explained in the text. The collective is not always 
more reversible as shown by the blue circles that correspond 
to a collective machine with colder baths. 


The last term scales like I/A^^, so the leading term is 
the coherence difference that scales like 1 /N. Since up 
to level permutation pdnr,i is equal to pw,N-\-i-i h follows 
that C{pdnr,i) = C{pw,N-\-i-i) and therefore 


^^bathsd T ^^baths,N-\-l—'i 


^{Pw,i\\Peq,i^ 

■T D(p^^ 7V'-|-1—i I |Peg, V+l—i) 

= 0{1/N^). (24) 


This pairwise cancellation of the leading order, I/A^, in 
the bath entropy generation, enables the reversible limit 
of the O = TT machine. After providing the mathematical 
explanation for this pairwise cancellation, it instructive 
to examine the unique physical circumstances that en¬ 
ables this result. The pairwise cancellation is a manifes¬ 
tation of a very interesting symbiotic relation between the 
engines in the upper hemisphere and the engines in the 
lower hemisphere. The upper units and lower units oper¬ 
ate in a very different manner. The upper ones produce 
work and coherence. The engines in the lower hemisphere 
also produce work but since they do not have population 
inversion, they must use coherence in order to produce 
work (see Eig. 3b). The coherence they need comes from 
the engines in the upper hemisphere. This shows what 
the lower hemisphere units gain from the upper ones. 
What about the other direction? Do the upper ones get 
something in return from the lower ones? The answer 
is yes. As explained earlier, and as shown in the graph 
in Eig. 7), due to the coherence injection the engines in 
the lower hemisphere have negative bath entropy genera¬ 
tion. Equation (24) implies that for large N this negative 


contribution compensates the positive bath entropy gen¬ 
eration of the upper hemisphere engine. In summary, 
the upper units ’sponsor’ the lower units, while the lower 
units clean up after the upper ones (see lower table in 
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Fig. 7). To the best of our knowledge this is the first ex¬ 
ample of symbiotic relations between different ’species’ 
of machines in thermodynamics. 

If N is large enough this machine operates as a re¬ 
frigerator when the order of the swap operations is re¬ 
versed. With some modifications that will be described 
elsewhere, the reverse operation can be improved. 

Note that there is no Carnot theorem for multiple 
baths. There are efficiency limits but for the same mul¬ 
tiple baths there are different reversible machines that 
operate at different efficiencies. Only in the two bath 
case the reversibility constraint is sufficient to fix the ef¬ 
ficiency. Hence, the ft = n heat engine has a reversible 
limit with effieieney rj = 1 — AEc/AE^ that is not fixed 
by the baths temperatures. 

VIII. CONCLUDING REMARKS 

In this work we have introduced two quantum ther¬ 
modynamics operations - coherence extraction (CE) and 
coherence injection (Cl). The first (CE) enables to reduce 
the entropy pollution while the latter can boost perfor¬ 
mance and extend the parameter regime in which the 
machine operates as an engine. By allowing heat ma¬ 
chines to interact with each other according to these two 
interactions, we studied what is, to the best our knowl¬ 
edge, the first example of coherent collective heat ma¬ 
chine. Eor certain parameters the work in this device 
scales quadratically with the number of composing en¬ 
gines and on top of the work boost it can also be more 
reversible compared to standalone operation. In addi¬ 
tion, the collective machines constructed here make use 
of coherence extraction and injection without external co¬ 
herence storage. The terminals of the collective machine 
are just thermal baths and a classical driving field (work 
repository). Nevertheless the internal coherence manipu¬ 
lation enable to extract more energy from the heat baths 
compared to individual standalone operation. 

The harvesting of coherence has a significant impact 
on the role of the second law in heat machine operation, 
as less information is erased in each cycle. This is neatly 
expressed in a simple relative entropy expression for the 


total entropy production. The relative entropy scaling 
shows that the information salvaged in the coherence can 
be significant enough to enable new reversible limits for 
otherwise irreversible machines. 

Eor a standard quantum machine, strong thermal- 
ization destroys coherence and turns the device into a 
stochastic thermal machine. The present study breaks 
this seemingly necessary connection between full ther- 
malization and stochastic-incoherent operation. In the 
collective operation mode, full thermalization and coher¬ 
ence effects coexist. 

The use of full swap for CE and Cl was chosen to avoid 
entanglement generation between particles as it may in¬ 
crease the entropy of the reduced systems and increase 
the entropy pollution. Nonetheless, it is possible that en¬ 
tanglement created by partial swap can be put to use in 
some variation of the engines studied here. 

CE and Cl give rise to new types of quantum thermo¬ 
dynamic devices. Experimental realizations of the ideas 
presented here call for techniques developed for quan¬ 
tum computation in order to implement the swap gates. 
In addition, open system techniques are needed to gener¬ 
ate the thermalization stroke. Although these techniques 
are not yet perfectly integrated together at the moment, 
both are feasible with present day technology. In par¬ 
ticular, it seems likely that the rapidly growing maturity 
of superconducting circuits technology will play a major 
role in quantum thermodynamics devices. In the context 
of the present setup, superconducting technology is well 
suited for fabricating sets of qubits that interact with 
their neighbors as well as with other baths [7]. 
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